Coherent pumping of a Mott insulator: Fermi golden rule versus Rabi oscillations 



0^ 
O 
O 
(N 

D 

\o 



I 

> 

e 

I 

o 
o 



> 

m 
in 
o 

oo 
O 



F. Hassleri and S. D. Huber^'i 

^Institute for Theoretical Physics, ETH Zurich, 8093 Zurich, Switzerland 
'^Department of Condensed Matter Physics, The Weizmann Institute of Science, Rehovot, 76100, Israel 

(Dated: February 26, 2009) 

Cold atoms provide a unique arena to study many-body systems far from equilibrium. Further- 
more, novel phases in cold atom systems are conveniently investigated by dynamical probes pushing 
the system out of equilibrium. Here, we discuss the pumping of doubly-occupied sites in a fermionic 
Mott insulator by a periodic modulation of the hopping amplitude. We show that deep in the 
insulating phase the many-body system can be mapped onto an effective two-level system which 
performs coherent Rabi oscillations due to the driving. Coupling of the two-level system to the 
remaining degrees of freedom renders the Rabi oscillations damped. We compare this scheme to an 
alternative description where the particles are incoherently pumped into a broad continuum. 

PACS numbers: 03.75.Ss, 71.10.Fd, 31.15.aq 



Cold atoms in optical lattices provided us with nu- 
merous insights into the paradigmatic quantum phase 
transition displayed by interacting lattice bosons [l], . 
Progress in the same setup with fermionic atoms has now 
led to new, equally-exciting results in the study of the 
Hubbard model I, H i, i, Q. A decrease of the compress- 
ibility when crossing over from the weakly to the strongly 
interacting region was reported [1,0. Furthermore, in 
the experiment by Jordens et al. [8[, the strongly cor- 
related regime of fermions in an optical lattice has been 
investigated by harmonically modulating the strength of 
the light beams forming the optical lattice. The modula- 
tion of the lattice generates doubly-occupied sites (dou- 
blons). This generation of doublons deep in the Mott 
insulator, where the interaction U dominates over the 
hopping to, is the focus of our present work. Experimen- 
tally, the fraction of particles on doubly-occupied sites 
(or double occupancy T)) was measured after modulat- 
ing the optical lattice. For U/6tQ = 13.6, a saturation 
was observed at around I?sat = 30% after a modulation 
time T„iod ~ 50h/U (50 "cycles") [8] which sets an upper 
bound on the saturation time Tsat < Tmod [see Fig.[TJa)]. 

We consider fermionic atoms, such as '^'^K, subject to 
an optical lattice in the framework of the Hubbard model 
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at half filling and with the interaction U dominating over 
the hopping amplitude tg. Here, the fermionic operators 
cjg. create an atom in one of two possible hyperfine states 
a =t,i at site i and the sum (i,j) runs over nearest 
neighbors. The ground state of the system consists of one 
particle per site with the hopping to suppressed due to the 
strong on-site repulsion U. As the temperature T in the 
experiment is larger than the super-exchange coupling 
J = 'Hq/U [sI, no significant quantum correlations over 
the bonds are present and the ground state of the system 
is largely composed of singly-occupied sites (T <C U) with 
a random spin- direct ion on each site. A harmonic mod- 
ulation of the lattice with frequency lu k, U /h transfers 
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FIG. 1: (color online) (a) Sketch of an experimental sequence 
for the generation of double occupancy by modulating the 
optical lattice during a time r. After an initial buildup of 
doubly-occupied sites with a rate Tup, the fraction of parti- 
cles residing in doubly-occupied sites saturates at a value Dsat 
after a time Taat . (b) Illustration of the effective Hilbert space 
HcS consisting of one bond embedded in an array of sites. 
The bath degrees of freedom consist of all doublon-hole con- 
figurations in which the pair resides on a different bond than 
where it was created. 



energy to the system and (possibly) promotes a particle 
across the bond, thereby creating a doubly-occupied site 
next to an empty site, a doublon-hole pair. 

Despite the fact that the Hubbard model represents a 
many-body problem in a huge Hilbert space, we adopt a 
simple description via modeling a single (typical) bond 
[Fig-HJb)], and taking the remaining bonds into account 
as an effective bath. Two different ideas lend themselves 
to modeling the process of pumping doublon-hole pairs 
[see Fig. [5]: (a) the direct excitation of the bond via 
incoherent transitions into the doublon-hole continuum 
with a rate given by the Fermi golden rule 
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here. Eg denotes the energy of the ground state, p is 
the density of states of the doublon-hole continuum, and 
(/[iJniodls) is the matrix element arising from the mod- 
ulation of the lattice. Taking the modulation strength 
to be given by 5t and estimating the density of states 
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arbitrary time 

FIG. 2: Comparison of two possible descriptions: (a) A rate 
equation description where an incoherent pump rate [rj^p''^] 
into a broad doublon-iiole continuum and subsequent decay 
['^downl incoherent stimulated emission [Fup''^] lead to a 
saturation, (b) Coherent Rabi oscillations [with frequency 
Qr] between the ground state and a single exited state where 
decay out of the effective Hilbert space is introduced via cou- 
pling to a bath [rin(out)]. 



as p oc 1/tQ, we obtain the scaling of the pump rate 
^n?od ^ ^f^/to as a function of the hopping amplitude 
to and the strength of the modulation St. (b) Coher- 
ent Rabi oscillations in the isolated bond between the 
ground state of the bond and the doubly-occupied ex- 
cited state with a Rabi frequency Qj^ oc 6t proportional 
to the modulation strength. The Rabi oscillations gen- 
erate doublon-hole pairs and lead to a periodic modula- 
tion (with frequency fJ^) of the double occupancy. Note 
that the rates of the two (extreme) processes, incoherent 
pumping into a band (a) (with rate T^"^) and coher- 
ent excitation of doublon-hole pairs (b) (with rate hil^f), 
have a different dependence on the modulation strength 
St and can therefore be discriminated experimentally. 

Here, we propose a third way where the coherent 
Rabi oscillations produce a doublon-hole pair and, sub- 
sequently, hopping of the doublon (or the hole) to the 
neighboring sites leads to an incoherent decay of the ex- 
citation into the doublon-hole continuum. The process is 
governed by two time scales, the Rabi frequency oc dt 
which periodically modulates the double occupancy and 
the decay Fhop esc to which leads to a saturation. 

A similar system, where a two-level atom is coupled 
simultaneously to a laser field and to vacuum modes of 
the radiation field, is well-studied in quantum optics p^ . 
The two-level atom corresponds to the two states on the 
bond and the laser plays the role of the harmonic mod- 
ulation of the lattice. However, in the atomic example, 
the spontaneously emitted photon takes the atom back 
into the ground state where it continues with the Rabi 
cycles, while in the present case the doublon (or hole) 
itself exits the system. The present process then resem- 
bles more strongly that of optical pumping into a third 
nondecaying level, with the accumulation in the pumped 
level playing the role of the buildup in double occupancy. 

In the following, we first describe the effective model in 
the single bond Hilbert space driven by a harmonic mod- 
ulation of the lattice and derive the frequency of the Rabi 
oscillations. We then account for the bath consisting of 
the surrounding sites via a master equation approach and 
analyze it in different limits. We compare this approach 
to the situation where the bond is directly excited into 



the doublon-hole continuum; see (a) above. It turns out 
that the time scale (r^) for the Rabi cycle is in rough 
agreement with the experimental value of the saturation 
time, whereas the treatment via a direct excitation into 
the doublon-hole continuum leads to a time scale (tfgr) 
which is an order of magnitude larger (r^^^ ^ Tj^ <^ 'T'fgr)- 
For a quantitative agreement, one has to go beyond a sin- 
gle bond in the description of the pumping process. Dif- 
ferent bonds are coupled via the pumping field, leading 
to a Dicke-like description [Tlj, but, most likely, many- 
body effects dominate the coupling of different bonds; a 
systematic study of larger clusters deserves future efforts. 

In deriving the effective Hamiltonian governing the 
dynamics in the Hilbert space of a single bond (the 
sites are labeled by L and R), we neglect the confin- 
ing potential and focus on a homogeneous three dimen- 
sional setup. The ground state at half-filling involves 
one particle per site without any spin order, i.e., the 
particles on the bond occupy the three triplet states 

and the singlet state | S* ) — -^[cl^cl^^—cl^ c]^^ ) 1 0) . These 
states constitute the four degenerate ground states of a 
single bond and are realized with equal probability of 1 /4 
each. A harmonic modulation of the lattice beam with 
V — Vq + 6V cos{u!t) leads to a time-dependent hopping 
Hamiltonian 
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with St /to = (3/4 - ^yVo/Er)SV/Vo and r denoting 
the time. This expression is valid for SV/Vq ^ 1 and 
SV/Vq <C \/Er/Vo 12]. For current experiments with 
fermions Q, where SV/Vq = 0.1 and Vq does not exceed 
lOEr, both conditions are safely fulfilled. In the deriva- 
tion of Eq. ([3|), we have dropped the constant hopping 
amplitude proportional to to, as hopping is blocked in the 
Mott insulator at half-filling. The hopping amplitude to 
will be reintroduced later as it is responsible for the de- 
cay of the doublon-hole pair by coupling the bond to the 
surrounding lattice. The modulation couples the ground 
states to the doublon-hole states |(iL(R)) = <^L(R)T'^L(R)il'^^' 
with the doublon located on the left (right) site. These 
states exhibit an energy offset U with respect to the en- 
ergy of the ground states. The total Hamiltonian in the 
Hilbert space spanned by the four ground states (triplet 
and singlet) and the two doublon-hole pairs is given by 

Has= U\da,)ido.\+H^op~ J2 U\d^){do.\ 

Q — L,R a — L,R 

+ St[e'^^/^\d,){S\ + e-'^^^^\S){d,\\, (4) 

where the first term describes the energy offset by U and 
the second term the driving due to the lattice modula- 
tion ([3]). In going from the first to the second equation 
in (III), the rotating wave approximation has been used 
and we have introduced the active and passive doublon- 



hole states given by \d. 



14) ± |dR))/V2; these are 
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related to the left (right) doublon-hole state via a basis 
change. Driving the system on resonance with huj = U , 
the singlet state is coupled to the excited state via Rabi 
oscillations. Within the present description, off-resonant 
coupling can be described in the standard manner. As 
can be seen from Eq. ((4]), the triplet states are unaffected 
by the lattice modulations. Only the singlet state is cou- 
pled to the active doublon-hole pairs |da). Going over to 
the interaction representation with the reference dynam- 
ics given by the first term in ([3]), the relevant Hamiltonian 
in the effective two-level system Heff = span{|da), 15*)} 
can be expressed as a matrix 



H 



5t 
St 



(5) 



which describes the coherent oscillations in the two-level 
system with the Rabi frequency 
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(6) 



The triplet states and the passive doublon-hole pair re- 
main cigenstates of the time-dependent Hamiltonian and 
therefore particles in such states are not driven. 

So far, we have neglected the effect of to on the excited 
state I da)- Provided the bond is in a doublon-hole state, 
we can lose the state from our single bond description 
via the hopping of the doublon (or the hole) out of the 
bond; cf. Fig.[l](b) This effect is captured by a bath 
representing the states with the doublon and hole sep- 
arated from each other. Note that the bath states still 
contribute to the double occupancy (as there is still a 
doublon present) but do not participate in Rabi oscilla- 
tions (as the doublon and the hole are not located on the 
same bond). We denote the rate at which the doublon- 
hole bond decays into a separated doublon and a hole 
by Tout, and the reverse rate at which a doublon and a 
hole collide to form a doublon-hole pair by Fin. Taking 
the coupling to the bath into account and performing a 
standard analysis similar to leads to an incoherent 
term (besides the standard term proportional to [H,p]) 
in the evolution of the reduced density matrix p in the 
effective Hilbert space Tioff — span{|(ia), \S)} of the form 



dp 
dr 



_ /^2(routPll — TinPo) routPl2 I (r,^ 

routP2i 



Different from a standard master equation [Tjl, Eq. ([7]) 
is not probability preserving, since with probability po = 
1 — Pii — P22 the doublon and the hole have separated 
and thus have left the effective Hilbert space HcS- In the 
steady state situation with drP = 0, the probability that 
the system contributes to the double occupancy is given 
by 
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(8) 



the sum of the probabilities that the system is in the 
doublon-hole pair state (pn) or in the bath (po). Note 
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FIG. 3: Fraction of particles residing in doubly-occupied sites 
after modulation time r. After a buildup time corresponding 
to approximately one 7r-pulse [tr = tt/JIr — h-K/25t], the 
decay, with (weak) strength Tout ~ 5t/h ~ 0.24to/fi and Fin — 
0, leads to a saturation at 'Dsa.t = 58% corresponding to the 
fact that all "active" bonds, singlets, are excited in the system 
(solid line). A more realistic situation where the damping 
is increased by a factor of 2 renders the oscillations barely 
visible (dashed line). Taking into account a finite value of 
Tin leads to a decrease of the saturation value. Note that, 
irrespective of the saturation value, the time scale needed to 
pump the system into a doubly occupied state is given by the 
Rabi frequency. 



that atoms which left the effective Hilbert space still con- 
tribute to the double occupancy. The bath can be pic- 
tured as a gas of doublons and holes (with equal num- 
bers). Initially, the bath does not contain any particles 
such that Fin = (this approximation remains valid as 
long as the density of particle is low); furthermore, we 
estimate the escape rate by Fout = Fhop ~ to/h. For 
Fin = 0, Eq. (HI) provides the saturation value p^^ — 1, 
i.e., the singlet is certainly transferred to one of the 
doublon-hole states. Using the dilute gas approximation, 
Fin = 0, we can solve the equation of motion for the 
density matrix 
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(9) 



Assuming that initially the bond is in a singlet state 

^0 



p{r = 0) = 



1 



(10) 



Eq. © gives the time evolution of the density matrix of 
the bond. Qualitatively, the weight is transferred to the 
exited state \da) at a rate Mlj^ and then the decay into 
other states occurs with a rate Fout (cf . Fig. [3]) . 

To calculate the fraction of particles in a doublon state, 
we need to know how many bonds are initially in a singlet 
configuration. We estimate the probability that a specific 
site is part of asinglet as Ps = 1 — (3/4)'* « 0.58, here d — 
3 is the "dimension" of the lattice modulation (it counts 
the number of bonds per site) and we have assumed that 
the singlets are independent of each other and occupied 
with a probability 1/4. The fraction of particles in a 
doublon state is given by 



T> = PsPu{t) 



(11) 



with Po {t) = 1 ~ P22 {t) the probability that this singlet 
is transferred to a doublon state. The system quickly 
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reaches a steady state with Vsat = Ps (cf- Fig. [S]). In 
three dimensions, we reach a maximal value I'sat — 58%. 
Note that at this value of double occupancy the dilute gas 
approximation Tm = becomes invalid. We expect that 
Fin increases for increasing double occupancy and, there- 
fore, the double occupancy saturates at a lower value. 
A theory that aims at a quantitative result for the sat- 
uration value should take such effects into account and 
solve the system self-consistently. Assuming that for long 
times the density of particles in the bath is high such 
that Fin » Fout, we obtain the result that p^^* — 50% 
[cf. Eq. ©land Pgat = 29%, comparable to experimen- 
tal results [Sj]. The time scale needed to excite the sys- 
tem is approximately given by the time Tr = Tr/f^R of 
a TT pulse, irrespective of the saturation value obtained. 
Note that the finite line-width of the perturbing laser, 
the scattering between doublons, and an inhomogeneous 
environment for a single bond due to the trap are other 
processes which lead to a dephasing of the coherent single 
bond dynamics. As the effect of tg on the excited state 
leads already to an over-damped dynamics, these other 
effects are neglected here. 

Another approach to describe the saturation in the 
double occupancy is based on a rate equation for the 
doublon fraction. The rates at which excitations are gen- 
erated are taken into account via a rate equation for the 
doublon concentration 

a,i? = (i-i?)rup-OTdown. (12) 

Within this description two independent processes are 
assumed. First, the system is excited with a rate Fup. 
Second, at a later uncorrelated stage, these excitations 
are removed with a rate Fdown [cf. Fig.[2{a)]. In order to 
compare the two schemes, the rate equation and the Rabi 
oscillation, we extract the time scales of both processes. 
The Fermi golden rule matrix element is approximately 
given by bt^ @ . Additionally, we need an estimate for the 
density of states. Once excited, the doublon as well as 
the hole can move freely through the lattice. The corre- 



sponding doublon-hole continuum acquires a width of ap- 
proximately 24io (in three dimensions). A crude estimate 
for the density of states is then given by p « l/24to- We 
then arrive at the two time scales Tr = tt/Or = 7rS/2(5t 
and TpGR = (ft/27r)24to/(5i2 = (2\m ■k'^ 8t)T^ « IOtr, as- 
sociated with a Rabi oscillation and the rate equation, 
respectively. For the last step, we have used (St/^o ~ 0.24 
corresponding to current experiments. Note that the 
time scale to reach saturation in the incoherent Fermi 
golden rule approach is an order of magnitude larger 
than the one predicted by the Rabi oscillation picture. In 
Ref. P the saturation sets in before 50 "cycles", which 
is even smaller than Tr [TVnod 

= hQhjU = (200(5t/C/)rR K. 
O.Gtr for U/&ta = 13.6 and 5t/tQ « 0.24], promoting 
the picture of the coherent Rabi oscillations. In an ex- 
periment, one can test one theory against the other by 
changing the coupling strength 5t while keeping all other 
parameters, such as the bandwidth l/24to, constant. The 
time scale for the saturation scales as 5t~'^ in the Fermi 
golden rule picture, whereas Rabi oscillations produce a 
saturation time inversely proportional to 5t. 

In summary, we have discussed two alternative expla- 
nations for the buildup and the saturation of the dou- 
ble occupancy in the modulation experiments of Ref. Q 
deep in the Mott insulating phase. The comparison to 
experiment favors the description in terms of an effec- 
tive two-level system with coherent driving, leading to 
the characteristic time scale of Tr k, tt/JIr = fnr/2St for 
half a Rabi cycle. Given the large decay rate propor- 
tional to to > h^ji « St, only one 7r-pulse is observable. 
The predictions can be experimentally tested by studying 
the dependence of the saturation time on various system 
parameters as Vo/Er, SV/Vq, and to/U. 
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